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13 Many hydrological models have been introduced in the hydrological literature to predict runoff 

14 but few of these have become common planning or decision-making tools, either because the data 

15 requirements are substantial or because the modeling processes are too complicated for operational 

16 application. On the other hand, progress in regional or global rainfall-runoff simulation has been 

17 constrained by the difficulty of measuring spatiotemporal variability of the primary causative factor, i.e. 

18 rainfall fluxes, continuously over space and time. Building on progress in remote sensing technology, 

19 researchers have improved the accuracy, coverage, and resolution of rainfall estimates by combining 

20 imagery from infrared, passive microwave, and space-borne radar sensors. 

21 Motivated by the recent increasing availability of global remote sensing data for estimating 

22 precipitation and describing land surface characteristics, this note reports a ballpark assessment of 

23 quasi-global runoff computed by incorporating satellite rainfall data and other remote sensing products 

24 in a relatively simple rainfall-runoff simulation approach: the Natural Resources Conservation Service 

25 (NRCS) runoff Curve Number (CN) method. Using an Antecedent Precipitation Index (API) as a 

26 proxy of antecedent moisture conditions, this note estimates time-varying NRCS-CN values 

27 determined by the 5-day normalized API. Driven by multi-year (1998-2006) Tropical Rainfall 

28 Measuring Mission (TRMM) Multi-satellite Precipitation Analysis, quasi-global runoff was 

29 retrospectively simulated with the NRCS-CN method and compared to Global Runoff Data Centre data 

30 at global and catchment scales. Results demonstrated the potential for using this simple method when 

31 diagnosing runoff values from satellite rainfall for the globe and for medium to large river basins. This 

32 work was done with the simple NRCS-CN method as a first-cut approach to understanding the 

33 challenges that lie ahead in advancing the satellite-based inference of global runoff. We expect that the 

34 successes and limitations revealed in this study will lay the basis for applying more advanced methods 

35 to capture the dynamic variability of the global hydrologic process for global runoff monitoring in real 

36 time. The essential ingredient in this work is the use of global satellite-based rainfall estimation. 
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Abstract 


2 Motivated by the recent increasing availability of global remote sensing data for estimating 

3 precipitation and describing land surface characteristics, this note reports a ballpark assessment of 

4 quasi-global runoff computed by incorporating satellite rainfall data and other remote sensing products 

5 in a relatively simple rainfall-runoff simulation approach: the Natural Resources Conservation Service 

6 (NRCS) runoff Curve Number (CN) method. Using an Antecedent Precipitation Index (API) as a 

7 proxy of antecedent moisture conditions, this note estimates time-varying NRCS-CN values 

8 determined by the 5-day normalized API. Driven by multi-year (1998-2006) Tropical Rainfall 

9 Measuring Mission (TRMM) Multi-satellite Precipitation Analysis, quasi-global runoff was 

10 retrospectively simulated with the NRCS-CN method and compared to Global Runoff Data Centre data 

1 1 at global and catchment scales. Results demonstrated the potential for using this simple method when 

12 diagnosing runoff values from satellite rainfall for the globe and for medium to large river basins. This 

13 work was done with the simple NRCS-CN method as a first-cut approach to understanding the 

14 challenges that lie ahead in advancing the satellite-based inference of global runoff. We expect that the 

15 successes and limitations revealed in this study will lay the basis for applying more advanced methods 

16 to capture the dynamic variability of the global hydrologic process for global runoff monitoring in real 

17 time. The essential ingredient in this work is the use of global satellite-based rainfall estimation. 

18 Key words: Rainfall-Runoff Modeling, Remote Sensing Precipitation, TRMM 
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1 1. Introduction 


2 Many hydrological models have been introduced in the hydrological literature to predict runoff 

3 (Singh, 1995) but few of these have become common planning or decision-making tools (Choi et al., 

4 2002), either because the data requirements are substantial or because the modeling processes are too 

5 complicated for operational application. On the other hand, progress in regional or global rainfall- 

6 runoff simulation has been constrained by the difficulty of measuring spatiotemporal variability of the 

7 primary causative factor, i.e. rainfall fluxes, continuously over space and time. Building on progress in 

8 remote sensing technology, researchers have improved the accuracy, coverage, and resolution of 

9 rainfall estimates by combining imagery from infrared, passive microwave, and space-borne radar 

10 sensors (Adler et al., 2003). Today, remote sensing imagery acquired and processed in real time can 

1 1 provide near-real-time rainfall at hydrologically relevant spatiotemporal scales (tens of kilometers and 

12 sub-daily; Hong et al., 2005; Huffman et al., 2007; Joyce et al., 2004; Sorooshian et al., 2000; Turk and 

13 Miller 2005). Over much of the globe, remote sensing precipitation estimates is the only available 

14 source of rainfall information, particularly in real time. Correspondingly, remote sensing has 

15 increasingly become a viable data source to augment the conventional hydrological rainfall-runoff 

1 6 simulation, especially for inaccessible regions or complex terrains, because remotely sensed imageries 

17 are able to monitor precipitation and identify land surface characteristics such as topography, stream 

18 network, land cover, vegetation etc. Artan et al. (2007) demonstrated the improved performance of 

19 remotely sensed precipitation data in hydrologic modeling when the hydrologic model was re- 

20 calibrated with satellite data rather than gauge rainfall over four subbasins of the Nile and Mekong 

21 Rivers. 

22 Motivated by the recent increasing availability of global remote sensing data for estimating 

23 precipitation and describing land surface characteristics, this note attempts to obtain a ballpark 

24 assessment of global runoff by incorporating satellite rainfall data and other remote sensing products 

25 through a relative simple rainfall-runoff simulation approach: the United States Natural Resources 

o 



1 Conservation Service (NRCS) runoff Curve Number (CN) method (USDA, 1986; Burges et ah, 1998). 

2 Its simplicity is especially critical for the vast un-gauged regions and geopolitically trans-boundary 

3 basins of the world. Our effort is a first approach to understanding a challenging problem that lies 

4 ahead in advancing satellite-based global runoff monitoring. Thus, the use of NRCS-CN should not be 

5 construed as a call for replacement of other more advanced methods for rainfall-runoff simulation. We 

6 expect that the successes and limitations revealed in this study will lay the basis for applying more 

7 advanced methods to capture the dynamic variability of the hydrologic process for global runoff 

8 monitoring in real time. The essential ingredient in this work is the use of global satellite-based rainfall 

9 estimation. 

10 Although Ponce and Hawkins (1996) indicated that the NRCS-CN method is widely used in the 

1 1 USA and other countries, they also criticized it as a simplistic methodology to simulate the 

12 sophisticated hydrological system. As an example, this method is imprecise for the monsoon type 

13 climate in Ethiopia (Mohammed et al., 2004). Taylor et al. (2006) also show that the annual runoff in 

14 the Volta river basin is a linear function of cumulative rainfall during the wet season when more than 

15 approximately 700 mm of rain has fallen. In a literature review, Choi et al. (2002) concluded that 

16 NRCS-CN has useful skill because it responds to major runoff-generating properties including soil 

17 type, land/use/treatment, and soil moisture conditions. They point out that it has been successfully 

18 applied to situations that include simple runoff calculation (Heaney et al., 2001), assessment of long- 

19 term hydrological impact on land use change (Harbor, 1994) for tens of years, stream-flow estimation 

20 for watersheds with no stream flow records (Bhaduri et al., 2000), and comprehensive 

21 hydrologic/water quality simulation (Srinivasan and Arnold, 1994; Engel, 1997; Burges et al., 1998; 

22 Rietz and Hawkins, 2000). Recently, Curtis et al. (2007) used satellite remote sensing rainfall and 

23 gauged runoff data to estimate CN for basins in eastern North Carolina. Harris and Hossain (2007) 

24 found that simpler approaches such as the NRCS-CN method to be more robust than more complicated 

25 schemes for the levels of uncertainty that exist in current satellite rainfall data products. On the other 
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1 hand, we note the risks of implementing this, or any other method without fully understanding its 

2 associated ‘uncertainty’. As such, we adopt the NRCS-CN method to estimate a first-cut global runoff 

3 by taking advantage of the first 9 years of rainfall estimates from the Tropical Rainfall Measuring 

4 Mission (TRMM) Multi-satellite Precipitation Analysis (TMPA; Huffman et al., 2007). 

5 In this note we first develop spatially distributed and time-variant CN maps for the global land 

6 surface. Driven by multi-year remote sensing rainfall, the NRCS-CN method is then used to compute 

7 the surface runoff for each grid independently and subsequently route the surface runoff to the 

8 watershed outlet through downstream cells (USACE, 2000). Finally, simulated quasi-global runoff is 

9 evaluated with Global Runoff Data Center (GRDC) observed runoff (Fekete et al. 2000) and Water- 

10 Balance-Model simulated runoff (Thomthwaite and Mather 1955; Steenhuis et al., 1986 ;_V o ro sm arty 

11 etal. 1998). 

12 

13 2. Mapping NRCS-CN 

14 2.1 Data 

15 The data sets (i.e., precipitation, soil information, and land cover) required by the NRCS-CN 

16 runoff generation scheme are all available globally with a well-established record in Earth system 

17 analysis (Fekete et al., 2000). Information on soil properties is obtained from the Digital Soil of the 

18 World published in 2003 by Food and Agriculture Organization of the United Nations 

19 (http : II www . fao . or g/ AG/ agl/ agll/ dsmw , him) . The Moderate Resolution Imaging Spectroradiometer 

20 (MODIS) land classification map is used as a surrogate for land use/cover, with 17 classes of land 

21 cover according to the International Geosphere-Biosphere Programme classification (Fridel et al., 

22 2002). Routing information is taken from the HYDRO lk ( http://lpdaac.usgs . gov/gtopo3 0/h vdro/) , 

23 which provides global coverage of topography such as elevation, slope, and flow direction etc. These 

24 geo-referenced datasets are of value for users who need to run hydrologic models on both regional and 

25 global scales. The rainfall data used in this study are from the NASA TMPA (Huffman et al., 2007; 

5 
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http://trmm.gsfc.nasa.gov) and the runoff data are from GRDC/University of New Hemisphere 
(http://www.grdc.sr.unh.edu/) . 


3 

4 2.2 Mapping NRCS-CN 

5 The NRCS-CN estimates surface runoff as a function of precipitation, soil type, land cover, and 

6 antecedent moisture conditions. The latter three factors are usually approximated by one parameter, the 

7 CN (USDA, 1986). In this case, the set of Equations (1-2) is used to partition rainfall into runoff and 

8 infiltration. 


9 

0= 

(P -IA + PR) 

(1) 

10 

PrP 5 - 400 - 254 

CN 

(2) 


1 1 where P is rainfall accumulation (mm/day); IA is initial abstraction; O is runoff generated by P; PR is 

12 potential retention; CN is the runoff curve number, with higher CN associated with higher runoff 

13 potential; and IA was approximated by 0.2 PR. 

14 CN values are approximated from the area's hydrologic soil group (HSG), land use/cover, and 

15 hydrologic condition, the two former factors being of greatest importance in determining its value 

16 (USDA, 1986). First, following the USDA (1986) handbook, a global HSG map is derived from the 

17 digital soil classification which includes 13 textural classes, an important indicator for infiltration rate 

18 (Table 1). Given “fair” moisture condition (defined below), the MODIS land cover classification and 

19 the HSG map are used to estimate CN by indexing into the standard lookup tables in USDA (1986) and 

20 NEH-4 (1997). Figure 1 shows the estimated climatological global CN map for fair moisture 

21 conditions, with higher value associated with larger mnoff potential. Thus, for a watershed on a coarse 

22 grid, a composite CN can be calculated as: 
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(3) 


CN = 


HACn, 
' S4 


2 in which CN com is the composite CN used for runoff volume computations; / = the index of subgrids or 

3 watershed subdivisions. Ai = the drainage area of area i. The composite CN values for several 

4 watersheds are listed in Table 2. 


6 3. Time-variant NRCS-CN and Runoff Simulation 

7 3.1 Time-variant NRCS-CN 

8 Note that the CN values displayed in Fig. 1 are for the “fair” hydrologic condition from standard 

9 lookup tables, which are used primarily for design applications. However, for the same rainfall amount 

10 there will be more runoff under wet conditions than under dry. In practice, lower and upper enveloping 

1 1 curves can be computed to determine the range of CN according to the Antecedent Moisture 

12 Conditions (AMC): 

CN" 


13 


14 


20 


CN = ■ 


2.281 -0.01 28 \CN" 


CN" 1 = 


CN" 


(4) 


(5) 


0.427 + 0.00573CTV/ 7 

1 5 where upper subscripts indicates the AMC, 1 being dry, II normal (average), and III wet (Hawkins 

16 1993). The change of AMC is closely related to antecedent precipitation (NEH-4, 1997). We apply the 

17 concept of an Antecedent Precipitation Index (API) to provide guidance on how to estimate the 

1 8 variation of CN values under dry or wet antecedent precipitation conditions. Kohler and Linsley (1951) 

19 define API as: 


API^Y j P l k- t (6) 


/=-] 
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1 where T is the number of antecedent days, k is the decay constant, and P is the precipitation during day 

2 t. The model is also known as “retained rainfall” (Singh 1989). Decay constant k is the antilog of the 

3 slope on a semi-log plot of soil moisture and time (Heggen, 2001). API practice suggests that k is 

4 generally between 0.80 and 0.98 (Viessman and Lewis 1996). Here we use decay constant k as 0.85 for 

5 demonstration purpose. API generally includes moisture conditions for the previous five days (or 

6 pentad; NEH-4, 1997). In order to obtain time-variant CN, the site-specified API is first normalized as: 

7 NAPI =■== (7) 

Pf>- f 

(=- 1 

8 where T=5 for pentads, the numerator is API, and the denominator is a normalizing operator with two 

9 components: average daily precipitation P and the ]T/C f series. The “dry” condition is defined as 

10 NAPI < 0.33, the “wet” condition is defined as NAPI>3, and the intermediate range 0.33-3 is the 

1 1 “fair” hydrological condition. By definition, the surface moisture conditions are delineated as dry (or 

12 wet) if any pentad API is less than one third (or larger than three times) of the climatologically 

1 3 averaged pentad API, and fair conditions for all others. Summarizing, the CN can be converted to dry, 

14 fair, or wet condition using Equations 4-7 according to the moisture conditions approximated by the 

15 pentad NAPI. 

16 Using the multi-year (1998-2006) satellite-based precipitation dataset from NASA TRMM, the 9- 

17 year climatological pentad API is shown in Figure 2a. Thus, given any date, the pentad NAPI can be 

18 determined and thus CN can be updated with Equations (4-7). For example, on Aug. 25 th , 2005, the 

19 pentad rainfall accumulation, pentad NAPI, resulting hydrological conditions (dry, fair, or wet), and 

20 the updated CN on the same date are shown in Figures 2b-e, respectively. 

21 

22 3.2 Runoff Simulation 
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1 Using the concept of NAPI and the NRCS-CN method (NEH-4, 1997), the TRMM-simuIated 

2 runoff (TRMM-CN) can be calculated and compared with three sets of GRDC annual climatological 

3 runoff fields: observed (OBS), Water Balance Model (WBM)-simulated, and composite (CMP) from 

4 the OBS and WBM (Fekete et al, 2000). The WBM used the water-balance model of Thorthwaite and 

5 Mather (1955) with a modified potential evaporation scheme from Vorosmarty et al. (1998), driven by 

6 input monthly air temperature and precipitation from Legates and Willmott (1990ab). Note that the 

7 three GRDC runoff climatologies span the period 1950-1979 with incomplete data records, while the 

8 TRMM-CN runoff is simulated for 9 years (1998-2006) of satellite rainfall with complete 

9 spatiotemporal coverage. One assumption here is that the change of rainfall between the two time 

10 periods is small enough so that the resulted runoff climatology is spatially consistent. Table 3 shows 

1 1 that TRMM-CN runoff corresponds more closely with WBM, having a relatively high correlation and 

12 low error. An intercomparison with the GRDC runoff observation demonstrates that the WBM has a 

13 moderate advantage over the TRMM-CN runoff: the correlation and root-mean-square-difference 

14 (rmsd) between GRDC OBS and WBM are 0.81 and 159.7 mm/year (or 0.437mm/day), respectively, 

15 which is slightly better than for the TRMM-CN case (Table 3). 

16 Figure 3a shows the annual mean runoff (mm/year) driven by TRMM daily precipitation for the 

17 same 9-year period, in comparison with the GRDC observed runoff climatology (Fig. 3b). Note that 

18 the gray areas indicate no data or water surface in figures 3a-b. By averaging areas covered by both 

19 TRMM-CN and GRDC runoff data, Figure 3c shows the TRMM-CN runoff zonal mean profile, 

20 against the OBS, WBM, and CMP. In general, the TRMM-CN zonal mean runoff more closely follows 

21 the three GRDC runoff profiles in the northern hemisphere than in the southern. We believe that this 

22 difference is the result of having many more samples in the northern hemisphere, as well as more- 

23 accurate GRDC data. Considering the TRMM-CN runoff difference as a function of basin area shows 

24 the TRMM-CN performance deviates more for basins smaller than 10,000 km 2 , with significantly 

25 better agreement for larger basins (Figure 4). 
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2 4 . Summary and Discussion 

3 Given the increasing availability of global geospatial data describing land surface characteristics, 

4 this note estimated a global CN map primarily based on soil property and land use/cover information 

5 under the “fair” moisture condition. Then using API as a proxy of AMC, this note further estimated 

6 time-variant CN values bounded by dry and wet AMC approximated by pentad normalized API. 

7 Finally, driven by satellite-based TMPA precipitation estimates, quasi-global runoff was simulated 

8 with the NRCS-CN method and was compared with GRDC runoff measurements for climatology and 

9 at the basin scale. 

10 Although we were able to demonstrate the potential for using the NRCS-CN runoff model when 

1 1 predicting ball-park runoff values from satellite rainfall for the globe and medium to large river basins, 

12 there remain several unanswered questions: First, among many methods to estimate CN values, 

13 Hawkins (1993) recognized that remote sensing data may not be adequate to define the “true” value of 

14 a CN. Thus, field surveys of basin characteristics should be conducted where feasible in order to obtain 

1 5 “true” soil and land cover data. Second, while this study recognized the uncertainty of the estimates of 

1 6 actual CN values and assumed that they likely fall within the enveloping wet (upper) and dry (lower) 

17 conditions approximated by the 5-day Normalized API, it may be possible to adjust the CN more 

18 precisely to account for local or regional information. Finally, one major unaddressed hydrological 

19 concern for rainfall-runoff applications of remotely sensed precipitation is the thorough evaluation of 

20 satellite-based rainfall estimation error and its nonlinear influence on rainfall-runoff modeling 

21 uncertainty in varying landscapes and climate regimes (Hong et al., 2006; Hossain and Anagnostou, 

22 2006; Villarini and Krajewski, 2007). Thus, while we conclude that this simple approach seems to 

23 provide a reliable tool when both the scale and estimation error of satellite data are large, we also urge 

24 similar studies using more sophisticated hydrological models, particularly seeking to serve the vast un- 

25 gauged regions and geopolitically trans-boundary basins of the world (Hossain et al., 2007). 
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Table 1 Hydrological Soil Group (HSG) derived from soil property 


HSG 

USD A SOIL 

TEXTURE 

CLASS 

SOIL CONTENTS 

% OF 

EARTH’S 

SURFACE 

PROPERTY 

A 

1,2,3 

sand, loamy sand or 
sandy loam types of 
soils 

4.69 

Low runoff potential and high infiltration 
rates even when thoroughly wetted; 
chiefly of deep, well to excessively 
drained sands or gravels 

B 

4,5,6 

silt loam, loam, or silt 

8.41 

moderate infiltration rate and consists 
soils chiefly with moderately fine to 
moderately coarse textures 

C 

7 

sandy clay loam 

3.98 

low infiltration rates when thoroughly 
wetted and consist chiefly of soils with 
moderately fine to fine structure. 

D 

8,9,10,11,12 

clay loam, silty clay 
loam, sandy clay, silty 
clay or clay 

5.78 

highest runoff potential, very low 
infiltration rates when thoroughly wetted 
and consist chiefly of clay soils 

0 

0 

Water bodies 

65.55 

N/A 

-1 

13 

Permanent ice/snow 

11.59 

N/A 





Table 2 The composite CN for several watersheds for “fair” hydrological conditions. 


Watershed 

Amazon 

Mississippi 

Yangtze 

Colorado 

Mekong 

Uruguay 

Sacramento 

Albany 

Composite CN 

75.484 

73.165 

81.787 

78.621 

62.355 

83.7 

77.425 

54.702 

Basin length(km) 

4327 

4184 

4734 

1807 

3977 

1424 

926 

951 

Area (km 2 ) 

5,853,804 

3,202,958 

1,794,242 

807,573 

773,737 

355,505 

192,563 

132,799 


Global surface-averaged CN=72.803 


Table 3 TRMM-CN runoff climatology in the latitude band 50°S-50°N) compared to GRDC observed 
(OBS), Water Balance Model (WBM), and the later two composite runoff (CMP). 



GRDC Runoff Climatology 

Statistics 

OBS 

WBM 

CMP 

Corr. Coef. 

0.75089 

0.79906 

0.78944 

Bias ratio 

1.2799 

1.1174 

1.1187 

Rmsd 

0.56mm/day 

0.48mm/day 

0.51mm/day 


















1 Figure Captions 

2 

3 Figure 1 . Global NRCS runoff Curve Number map derived from USDA Hydrological Soil Groups and 

4 Land Cover Classification for fair hydrological conditions. 

5 

6 Figure 2. (a) Climatological pentad Antecedent Precipitation Index (API) averaged over 9 years (1998- 

7 2006); (b) pentad antecedent rainfall accumulation (mm) ending on August 25, 2005; (c) pentad 

8 Normalized API (NAPI) on August 25, 2005; (d) hydrological condition, with -1,0, 1,2 

9 corresponding to no data, dry, fair, wet conditions, determined by NAPI as of August 25, 2005; 

10 and (e) the updated CN on August 25,. 2005. 

11 

12 Figure 3. (a) The annual mean runoff (mm/year) simulated using NRCS-CN methods from TMPA 

13 estimates for the period 1998-2006; (b) the GRDC observed runoff (mm/year); and (c) runoff 

14 zonal mean profiles comparing TMPA precipitation (green) and simulated runoff (red) to 

1 5 GRDC runoff (blue) from the observed (left), WBGS (center)_, and composite (right data sets.. 

16 Note: the gray areas in (a) and (b) indicate no data or water surface. 

17 

18 Figure 4. TRMM-CN runoff difference distribution (a) and root-mean-square difference (b) as a 

19 function of basin area. 
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Figure 1. Global NRCS runoff Curve Number 
map derived from USDA Hydrological Soil 
Groups and land cover/use classification for fair 
hydrological conditions. 
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Figure 2. (a) Climatological pentad Antecedent Precipitation Index 
(API)) averaged over 9-year (1998-2006); (b) pentad antecedent 
rainfall accumulation (mm) ending on August 25, 2005; (c) pentad 
Normalized API (NAPI) on Aug 25, 2005; (d) hydrological condition, 
with -1, 0, 1,2 corresponding to no data, dry, fair, wet conditions, 
determined by NAPI as of August 25, 2005; and (e) the updated CN on 
date August 25, 2005. 
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Figure 3. (a) The annual mean runoff (mm/year) simulated using 
NRCS-CN methods from TMPA estimates for the period 1998- 
2006; (b) the GRDC observed runoff (mm/year); and (c) runoff 
zonal mean profiles comparing TMPA precipitation (green) and 
simulated runoff (red) to GRDC runoff (blue) from the observed 
(left), WBGS (center), and composite (right data sets.. Note: the 
gray areas in (a) and (b) indicate no data or water surface. 
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Figure 4. TRMM-CN runoff difference distribution (a) and root-mean-square 
difference as a function of basin area. 


